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I.  INTRODUCTION 


The  use  of  carbon  materials  in  advanced  thermal  protection  system 

applications  requires  an  understanding  of  the  thermochemistry  of  carbon  at 

o  1  2 

temperatures  exceeding  4000oKA’  .  Of  primary  importance  is  an  accurate 
knowledge  of  the  melt  temperature,  the  species  vapor  pressures,  and  species 
vaporization  coefficients,  since  these  properties  directly  affect  the  energy- 
dissipating  efficiency.  At  present,  the  uncertainty  in  the  melt  temperature 
of  carbon  is  at  least  600°K,  and  recommended  thermochemical  data  from  dif¬ 
ferent  investigators  result  in  predicted  total  vapor  pressures  differing  by 

o  3-5 

more  than  an  order  of  magnitude  above  4000  K  .  Of  equal  importance  is 
an  understanding  of  the  nonequilibrium  nature  of  the  phase  change  process  when 
massive  interphase  mass  transfer  takes  place  such  as  that  induced  by  high 
intensity  thermal  radiation. 


1  2 

The  approach  of  our  work  ’  has  been  to  address  simultaneously  the 
needs  of  high-temperature  carbon  thermochemical  property  data  and  nonequili¬ 
brium  phase  change  modeling.  The  development  of  a  nonequilibrium  phase  change 
model  first  allows  the  interpretation  of  laser-induced  phase  change  experi¬ 
mental  data;  i.e.,  equilibrium  thermochemical  property  information  can  be 

1  - 

"Proposal  to  Perform  a  Carbon  Nonequilibrium  Phase  Change  Research  Program," 
The  Aerospace  Corporation,  El  Segundo,  Calif.,  April  1979. 

2 

"Proposal  to  Continue  the  Carbon  Nonequilibrium  Phase  Change  Research  Pro¬ 
gram,"  The  Aerospace  Corporation,  El  Segundo,  Calif.,  August  1979. 

3 

Baker,  R.L.,  "An  Irreversible  Thermodynamics  Model  for  Graphite  Sublimation 
in  Radiation  Environments,"  Progress  in  Astronautics  and  Aeronautics:  Outer 
Planet  Heating  and  Thermal  Protection  Systems,  64,  R.  Viskanta  (ed.),  AIAA, 

New  York,  1575,  pp.  110-Z27: - - 

4 

Baker,  R.L.  and  P.G.  Crowell,  "Graphite  Material  Ablation  Performance  in  High 
Thermal  Radiation  Environments,"  Progress  in  Astronautics  and  Aeronautics: 
Entry  Heating  and  Thermal  Protection,  69,  W.  Olstad  (ed.),  AlAA,  New  Vork, 

I9SQ,  pp.  TSP22TT - 

5 

Baker,  R.L.,  "Graphite  Sublimation  Chemistry  Nonequilibrium  Effects,"  AIAA 
Journal,  15,  Oct.  1977,  pp.  1391-1397. 


obtained  from  a  nonequilibrium  experiment.  After  the  required  information  has 
been  obtained,  the  same  model  then  can  be  utilized  for  design  applications 
involving  carbon  materials  at  extreme  temperatures  and  under  nonequilibrium 
conditions. 

Toward  this  end  we  have  developed,  in  the  past  year,  a  nonlinear  non¬ 
equilibrium  model  applicable  to  the  carbon  phase  change  process  which  allows 

either  direct  sublimation  or  melting  with  subsequent  vaporization.  This  work 

3-5 

is  an  extension  of  our  earlier  work  involving  linearized  models  .  In 

addition,  we  have  obtained,  partly  in  collaborative  work  with  NASA  Ames,  high- 

temperature  (3500  to  4700°K)  carbon  ablation  data  using  laser  radiation  at 

2 

power  densities  ranging  from  10  to  50  kW/cm  in  continuous  wave  laser  ex- 

2  2 

periments  and  from  200  kW/cm  up  to  4  MW/cm  in  pulsed  laser  experiments. 

The  major  overall  objective  of  this  work  is  to  use  the  model  develop¬ 
ment  effort,  discussed  in  Sections  II  and  III,  to  interpret  the  laser-induced 
carbon  nonequilibrium  phase  change  data,  described  in  Sections  IV  and  V,  in 
such  a  way  as  to  obtain  the  needed  information,  i.e.,  carbon  species  heats  of 
formation,  free  energy  functions,  and  vaporization  coefficients,  as  well  as 
the  melt  temperature  of  carbon.  Preliminary  data  analyses  are  included  as 
part  of  the  continuous  wave  laser  and  pulsed  laser  data  discussions  in  Sec¬ 
tions  IV  and  V,  respectively.  The  melt  temperature  question  is  discussed  in 
Section  VI.  Detailed  analysis  of  these  data,  to  obtain  a  critical  evaluation 
of  the  validity  of  the  JANNAF  thermochemical  property  data  and  the  melt  tem¬ 
perature  of  carbon,  is  planned  as  part  of  the  continuation  of  this  work. 


II.  MELT  LAYER  MODEL 


The  surface  boundary  condition  for  the  heat  conduction  into  the  bulk 
material  becomes  fundamentally  different  when  a  semitransparent  melt  layer 
occurs  on  the  surface  and  the  energy  input  is  supplied  by  radiation.  Thus,  in 
order  to  compare  experimental  pulsed  laser  data  with  results  predicted  from 
the  nonlinear,  nonequilibrium  phase  change  analysis  described  in  Section  III, 
proper  accounting  of  this  melting  surface  boundary  condition  must  be  made.  An 
equivalent  melt  layer  problem,  where  the  energy  input  is  supplied  by  a  heat 
source  or  forced  convection  heat  transfer,  was  described  by  Landau®  and  more 
recently  by  Crowell^.  Their  approach  is  used  below  to  describe  the  case  in 
which  simultaneous  radiation  absorption  and  thermal  conduction  occur  in  the 
melt  layer.  The  problem  formulation  is  given  in  Subsection  I I - A ,  and  repre¬ 
sentative  calculated  results  and  applications  are  discussed  in  Subsection 
II-B.  The  model  must  be  used  in  conjunction  with  the  nonlinear  Knudsen  layer 
models  described  in  Section  III  for  a  1-1  data  analyses  in  which  carbon  melting 
at  the  surface  is  postulated  to  have  occurred. 

A.  PROBLEM  FORMULATION 


1.  Temperature  Distributions 


The  problem  considered  is  shown  schematically  in  Fig.  1.  Laser  radia¬ 
tion  of  intensity  IQ  impinges  upon  the  surface.  It  is  assumed  that  IQ  is 
sufficiently  high  such  that  melting  of  the  material  occurs,  with  a  melt  layer 
of  thickness  6  formed  at  the  surface  above  the  melt  isotherm.  Within  the 
melt  layer,  laser  radiation  energy  is  simultaneously  absorbed  and  transported 
by  thermal  conduction  due  to  temperature  gradients.  Following  Refs.  6.  and  7, 
the  ablation-melting  process  is  assumed  to  be  steady  state; 

l 

Landau,  H.G.,  "Heat  Conduction  in  a  Melting  Solid,"  Quart.  Appl.  Phys.  36, 
Feb.  1965,  pp.  462-468. 

7 

Crowell,  P.G.,  "The  Nonequilibrium  Ablation  of  Carbon,"  Report  No.  TR-0079 
(4550- 76 ) -1 ,  The  Aerospace  Corporation,  El  Segundo,  Calif.,  Nov.  1978. 


Steady  State  Melt  Layer 


therefore,  the  melt  layer  surface  and  the  melt  isotherm  remain  separated  by 
the  distance  0,  and  both  surfaces  translate  into  the  bulk  material  at  the 
steady-state  ablation  rate  s,  i.e.,  the  surface  velocity. 


Defining  one-dimensional  stationary  and  translating  coordinate  systems 
x'  and  x  where  x1  =  x  +  st,  we  find  the  energy  conservation  equation  for  the 
melt  layer  in  the  stationary  x'  coordinate  system  to  be 


dJi 


-  K I 


e-*(x'-st)=  k 


L 


(II-l) 


Following  Landau,  this  equation  may  be  transformed  into  the  moving  reference 
frame  with  coordinate  x  to  obtain 


•  dTL 
dx2  dX 


"LCP.  ^ 


-  kx 


Pl"l 


(II-2) 


where  <?^  is  the  thermal  diffusivity  of  the  liquid  and  k  is  the  spectral 
absorption  coefficient.  A  time-dependent  equation  equivalent  to  Eq.  ( 1 1 -2 ) 

o 

for  the  nonmelting  case  has  been  reported  by  Dabby  and  Paek  . 


In  the  present  analysis,  it  has  been  assumed  that  all  radiation  absorp 
tion  within  the  solid  takes  place  at  the  solid  surface.  Thus,  following  the 
procedure  described  above,  the  differential  equation  describing  the  solid  tern 
perature  in  the  region  below  the  melt  isotherm  is 


dTf 


s  J  _  n 

c:sr  - 0 


( 1 1-3 ) 


S - 

Dabby,  F.W.  and  Paek,  U.-C.,  "High-intensity  Laser  Induced  Vaporization  and 
Explosion  of  Solid  Material,"  IEEE  J.  of  Quantum  Electronics,  QE-8(2).  Feb. 
1972,  pp.  106-111. 
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The  boundary  conditions  for  these  second-order,  ordinary  differential  equa¬ 
tions  are 


At  x  *  0;  T.  =  T 
L  s 

At  X  =  6i  T,  =  T,  T<-  =  T 
L  m  5  m 

At  x  *  oo;  Ts  =  Tb  ( 1 1-4 ) 


where  T  ,  Tm,  and  are,  respectively,  the  surface  temperature  of  the 

melt  layer,  the  melt  temperature,  and  the  bulk  temperature  of  the  solid  (see 

Fig.  1).  In  all  cases  considered,  it  is  assumed  that  T  >  T  »  T.  . 

s  m  b 

By  direct  application  of  classical  linear  equation  solution  methods  and 
the  given  boundary  conditions,  the  following  expressions  are  obtained  for  the 
liquid  and  solid  temperature  distributions: 

t  lTe  -  TJ  e(‘s/aL)x  +  T  -  T  e('s/<rL)3 
rft  Is  ml  ms 

L  - ; — (wnr; - 

l>-»  I 

>o  r  1  1  |[l-e-|fi][l-e<4/‘,l’>]  ,  x  S  « 

pjTTT  [a'-IJ  -s  7oT&  | 

L  PL  L  J  1  -  e  L  ( II-5 ) 


T 


S 


e-s(x-6)/<r$  + 


x  >  6 


(II-6) 
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The  dimensionless  group  a'  appearing  in  £ q.  ( I I -5 )  is  given  by 


a '  U,<7l,s) 


#C<T|_ 

~r 


( 1 1-7  a) 


2.  Determination  of  Melt  Layer  Thickness  3 

While  Eqs.  ( 1 1 -6 )  and  ( 1 1 -7 )  give  the  liquid  layer  and  solid  temper¬ 
ature  distributions,  they  are  expressed  in  terms  of  the  unknown  thickness  of 
the  melt  layer  6.  To  evaluate  6,  it  is  necessary  to  supply  additional  infor¬ 
mation  which  may  be  obtained  from  the  energy  balance  equation  at  the  liquid- 
solid  interface.  In  obtaining  Eqs.  ( 1 1 -5 )  and  ( 1 1-6 ) ,  it  has  been  assumed 
that  material  property  values  are  constant.  In  the  following,  it  is  further 

assumed  that  p,  =  pc  =  p  and  C  =  C  =  C  but  not  that  k.  =  kc. 

L  5  PL  Ps  p  L  5 

The  energy  balance  equation  at  the  liquid-solid  interface,  x  *  6, 
assuming  unit  absorptivity  of  the  solid  surface  and  neglecting  emitted  energy 
at  temperature  T  ,  is 


Ioe-*-  qL(«)  -  qs(«)  =  ps  AHp 


( 1 1 -8 ) 


In  this  equation,  the  first  term  represents  the  radiation  energy  trans¬ 
mitted  through  the  liquid  layer  and  absorbed  at  the  surface  of  the  solid. 

Heat  conduction  due  to  temperature  gradients  is  given  by  qL  and  q^  for  the 
liquid  and  the  solid,  respectively.  The  balance  of  these  terms  supplies  the 
energy  to  melt  the  solid  at  the  interface  with  the  heat  of  fusion  given  by 
AHp.  The  terms  qL  (6)  and  q^  (6)  are  obtained  from  the  liquid  and  solid 
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temperature  distribution  Eqs.  ( 1 1-5 )  and  ( 1 1-6 ) .  Thus,  differentiating  these 
expressions,  we  obtain 
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(II-9) 


qs(«>  -  -Tb] 


(II- 10 ) 


Substituting  these  expressions  into  Eq.  ( 1 1-8 ) ,  we  obtain  the  following 
implicit  equation  for  the  evaluation  of  6 


A  +  BZ  +  CZ  °*  ’  -  0 


where  a' (*,aL,  s)  is  given  by  Eq.  (II-7a)  and 


(II-U) 


1 


B  =  ^  [Cp(Ts  -  V  *  AHf]  *  [?4t]  (H-12,i) 


c  = 


” 


'1  1 


(II-12c) 


z  *  J^*06 


( 1 1 -12  d ) 


For  specified  input  parameters,  Eq.  (.11-11)  may  be  solved  iteratively  for  Z. 
Then  6  is  obtained  directly  fr:m 


<5 


(II -12e) 


The  final  quantity  of  major  interest  and  the  main  reason  for  carrying 
out  the  present  analysis  is  the  value  of  the  heat  conduction  at  the  surface  of 
the  liquid  layer.  From  Eq.  (II-5)  we  obtain  by  differentiation 


dTL 

=  '  kL  Hx~ 
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(11-13) 
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As  discussed  in  the  calculated  results  in  Subsection  I I - B ,  this  term  is  always 
negative  so  that  radiation  energy  absorbed  in  the  liquid  layer  is  conducted 
back  to  the  liquid  surface  where  it  supplies  the  energy  required  to  vaporize 
the  liquid.  When  the  various  energy  transport  terms  are  utilized  as  shown 
schematically  in  Fig.  2,  it  can  be  shown  that  in  addition  to  the  liquid-solid 
interface  energy  balance  equation,  ( 1 1-8 ) ,  the  analytic  equations  given  here 
also  satisfy  the  overall  liquid  layer  energy  balance  given  by 


I0(l  -  e'K<5)  +  qL(0)  -  qL((5)  -  pi( Ts  -  Tm)  =  0  (11-14) 

the  global  energy  balance  equation 


lo  ‘  +  ^  AHF  +  %iv  =  0  (H-15) 

and  the  liquid  surface  energy  balance 


where 


(11-16) 


qT  *PSCP(TS  -  Tb),  q  AHp  =ps  AHf 


>S  L 


ps  AHy  +  C  (T, 


Ts) 


(H-17) 
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B.  CALCULATED  RESULTS  AND  APPLICATIONS 


1.  Temperature  Profiles 


The  temperature  T$  in  the  solid  below  the  liquid-solid  interface  is 
given  by  Eq.  ( 1 1 -6 ) ;  T^  decreases  monotonically  and  exponentially  from  the 
melt  temperature  Tm  at  the  interface  to  the  bulk  solid  temperature  Tp  far 
away  from  the  interface.  In  contrast,  the  liquid  melt  layer  temperature  T^, 
as  given  by  Eq.  (II-5),  exhibits  a  maximum  value  within  the  melt  layer  for  all 
positive  values  of  a' . 

Because  of  this  maximum  in  the  liquid  temperature  profile,  the  tempera¬ 
ture  gradient  dT^/dx  is  zero  where  T^  -  TL  ,  positive  at  x  =  0,  and 
negative  at  x  *  6 .  Thus,  radiation  energym£$sorbed  within  the  liquid  layer  is 
conducted  in  two  directions  away  from  the  plane  within  the  liquid  layer  having 
temperature  T^  .  At  the  surface  x  =  0,  the  heat  conduction  qL  (0) 

(Eq.  (11-16))  §8&pl  ies  the  energy  required  to  vaporize  the  liquid.  At  the 
surface  x  =  6,  the  heat  conduction  qL(d)  (Eq.  ( I I -9 ) )  supplies  a  portion  of 
the  energy  required  to  melt  the  solid  (see  Eq.  { 1 1-8 ) ) . 

For  specified  values  of  C  ,  T  ,  T  and  the  effective  heat  of  ablation 

p  iTr  S 

Q*  =  IQ/(ps),  the  temperature  profile  T^  =  TL(x/£)  is  only  a  function  of  a1. 
Using  the  effective  heat  of  ablation  Q*,  one  may  write  a'  alternatively  as 


a' 


K<r,pQ 


( I I - 7b) 


Also,  the  term  I0/(PLCn^  in  (H-5)  becomes  Q*/Cp.  For  an  assumed  melt 
temperature  T^  of  4200  K,  surface  tenperature  T$  of  4400°K,  Cp  of  5/12  cal/g-°K, 
and  Q*  obtained  directly  from  nonequilibrium  calculations  as  described  in 
Section  III,  the  calculated  liquid  layer  temperature  profiles  are  shown  as  a 
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function  of  a1  in  Fig.  3.  For  the  lower  range  of  ot1 ,  where  the  predicted 
value  of  becomes  much  larger  than  T$,  it  is  unlikely  that  the  pre¬ 

dicted  temp??lture  profile  is  achieved.  More  likely,  in  such  cases,  is  a 
nucleate  boiling  phenomenon  in  the  region  around  =  T^  or  possibly 
even  an  explosive  vaporization  phenomenon  such  as  that  sQ^ested  by  Dabby  and 

Q  Q 

Paek  for  solids.  For  values  of  T$  other  than  4400  K,  a  similar  depen¬ 
dence  of  the  temperature  profile  on  a'  occurs. 

2.  Melt  Layer  Thickness 

While  the  temperature  profile  T^(x/0)  depends  only  upon  ot1 ,  the 
liquid  layer  thickness  6  depends  upon  both  a^/s  and  k.  However,  the  depen¬ 
dence  upon  a,J%  is  given  explicity  by  Eq.  (II-12e),  and  the  quantity  Z  de¬ 
pends  only  upon  a'  through  the  implicit  function  represented  by  Eq.  (11-11). 
Calculated  values  of  Z  using  Eq.  ( 1 1 -1 1 )  with  T$  =  4400°K  and  = 

4200°K  are  shown  in  Fig.  4  as  a  function  of  o' .  Calculated  results  for 
other  values  of  Ts  ranging  from  4250  to  4500°K  also  have  been  made,  and 
the  dependence  of  Z  on  T$  is  very  weak. 

From  Fig.  4  and  Eq.  (II-7a),  the  dependence  of  melt  layer  thickness  6 

on  spectral  absorptivity  k  for  given  a^/i  can  be  readily  obtained.  Such 

results  are  shown  in  Fig.  5.  For  given  tf^/s,  6  decreases  by  a  factor  of  30 

as  k  is  increased  two  orders  of  magnitude.  For  all  the  calculations  carried 

out,  a'  varied  from  0.055  to  96.4,  and  the  predicted  liquid  layer  thicknesses 

varied  by  a  factor  of  about  70.  Thermal  diffusivity  of  the  liquid  was  varied 

an  order  of  magnitude  from  the  nominal  solid  carbon  value  of  0.93  x  10"^ 

?  -2  2 

cm  /sec  downward  to  0.93  x  10  cm  /sec.  The  recession  rate  also  was 
varied  about  an  order  of  magnitude  from  0.95  cm/sec  to  8.6  cm/sec. 
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Melt  Layer 


By  using  an  effective  heat  of  ablation  Q*,  we  may  express  a'  as  a  func 
tion  of  IQ,  i.e.,  Eq.  ( I I -7b) .  Then,  substituting  s  =  Iq/(PQ*)  into  Eq. 

(II -12e ) ,  we  obtain  the  following  equation  for  6  as  a  function  of  IQ: 

pQ*o ;  . 

6  = - p1  /n  2(  Of  [  IQ]  )|  (11-18) 

o 

4  3  -1 

The  dependence  of  6  on  I  for  pQ*  =  5.36  x  10  0/cm  ,  k  -  984  cm  ,  and 
-2  2  0 

<tl  =  1.86  x  10  cm  /sec  is  shown  in  Fig.  6.  The  melt  layer  thick¬ 
ness  decreases  about  20  percent,  while  I  increases  an  order  of  magnitude. 

The  corresponding  surface  temperature  T$  varied  from  4250  to  5000°K. 
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III.  NONLINEAR  KNUDSEN  LAYER  MODELS 


A.  BACKGROUND  AND  PREVIOUS  WORK 


1.  Problem  Description 


A  basic  premise  of  our  work  on  this  contract  since  its  inception  has 
been  that  carbon  high-temperature  (3500  to  4500°K)  equilibrium  thermo¬ 
chemical  properties,  i.e.,  species  vapor  pressures  and  melt  temperature,  can 
be  determined  from  nonequilibrium  phase  change  laser  heating  experi- 
ments  ’  .  The  success  of  such  an  approach  is  dependent  upon  modeling  non¬ 
equilibrium  behavior  in  the  Knudsen  layer  in  order  to  interpret  laser  heating 
experiments  in  terms  of  equilibrium  properties. 

The  essential  details  of  the  problem  are  illustrated  in  Fig.  7.  The 
Knudsen  layer  is  a  thin  region  (the  order  of  a  few  mean  free  paths  in  thick¬ 
ness)  inmediately  above  the  subliming  or  vaporizing  carbon  surface,  where  non¬ 
equilibrium  gas  dynamic  and  chemical  rate  processes  occur  which  determine  the 
change  in  the  carbon  vapor  state  as  it  moves  away  from  the  ablating  surface. 
Since  macroscopic  temperature,  pressure,  and  species  concentration  measure¬ 
ments  can  only  be  made,  at  the  present  time,  at  the  outer  edge  of  the  Knudsen 
layer,  a  Knudsen  layer  model  provides  a  link  between  these  measurable  macro¬ 
scopic  variables  and  the  equilibrium  (saturation)  thermodynamic  properties  of 
the  carbon  vapor  at  the  solid  (or  liquid)  surface.  Consideration  of  the 
Knudsen  layer  becomes  of  first  order  importance  for  phase  change  problems  in 

which  the  vapor  pressure  of  the  vaporizing  or  subliming  material  becomes 

3-5 

greater  than  the  ambient  pressure  .  This  occurs  typically  in  laser- 
induced  phase  change  processes. 

The  basic  output  of  a  Knudsen  layer  analysis  provides  the  “jump"  condi¬ 
tions  across  the  layer  and  thus  relates  the  gas-dynamic  properties  at  the 
Knudsen  layer  edge  to  equilibrium  thermochemical  properties  of  the  solid  or 


Knudsen  Layer 


liquid.  These  jump  conditions  have  been  described  as  being  analogous  to  the 

q  m 

Rankine-Hugoniot  jump  conditions  for  a  normal  gas-dynamic  shock  wave  ’  . 

However,  some  important  distinctions  must  be  kept  in  mind  when  speaking  of 
such  an  analogy.  First,  a  shock  wave  is  a  compression  region,  and  a  Knudsen 
layer  is  an  expansion  region.  The  most  important  difference  is  that  a  shock 
wave  may  be  described  mathematically  by  continuum  hydrodynamic  equations, 
unlike  a  Knudsen  layer.  The  reason  for  this  is  that  the  velocity  distribution 
function  across  a  Knudsen  layer  is  strongly  nonequilibrium  or  non-Maxwell i an. 
In  fact,  the  Knudsen  layer  edge  is  defined  as  the  location  at  which  the  velo¬ 
city  distribution  function  becomes  near  equilibrium.  Thus,  the  irreversible 
phenomena  which  determine  the  entropy  production  across  a  shock  wave  and 
across  a  Knudsen  layer  are  somewhat  different. 

2.  Linearized  Models 


Most  Knudsen  layer  analyses  carried  out  by  early  investigators 
neglected  the  momentum  equation,  and  sometimes  also  the  energy  equation,  and 
retained  only  those  terms  which  were  linear  in  the  jump  conditions.  An  excel 
lent  sunmary  of  approaches  utilizing  linearized  equations  is  given  by 
Wiechert^. 


4 

Baker  adapted  one  of  the  linearized  models,  which  utilized  a  non- 

12 

equilibrium  thermodynamics  approach  ,  to  the  carbon  phase  change  system  by 
assuming  that  carbon  vapor  is  made  up  of  species  molecules  only.  More 

1 - 

Ytrehus,  T.,  "Theory  and  Experiments  on  Gas  Kinetics  in  Evaporation,"  Pro¬ 
gress  in  Astronautics  and  Aeronautics:  Rarefied  Gas  Dynamics,  51(11),  J.L. 
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II,"  J.  Phys.  C:  Solid  State  Phys..  9,  1976,  pp.  553-569. 
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recently,  this  work,  retaining  linearized  equations,  was  extended  to  include 
carbon  species  through  and  to  consider  the  bounding  cases  of  frozen 
and  equilibrium  chemistry  in  the  Knudsen  layer^.  Comparison  of  calculated 

c 

results  using  this  model  with  pulsed  laser  graphite  ablation  data  led  to 
the  tentative  conclusions  that  carbon  melting  may  have  occurred  at  about 

4250°K,  and  that  the  JANNAF  carbon  species  C,  thermochemical  proper- 

13  J 

ties  may  not  be  correct. 

3.  Nonlinear  Models 


Recently,  Knudsen  layer  models  retaining  nonlinear  terms  in  the  mass, 

momentum,  and  energy  conservation  equations,  which  determine  the  jumps  in 

g  in 

intensive  variables  across  the  layer,  have  been  reported  *  .  This  work 

14 

appears  to  be  based  upon  the  earlier  pioneering  work  of  Anisimov  .  Analy¬ 
tical  work  carried  out  during  the  past  contract  year  has  consisted  of  extend- 

4  5 

ing  our  previous  linearized  models  ’  by  retaining  all  nonlinear  terms  in 
the  balance  equations.  Initial  work  consisted  of  extending  the  Anisimov  equa¬ 
tions  to  include  arbitrary  values  of  the  vaporization  coefficient  a  and  the 
specific  heat  ratio  y.  Anisimov's  results  were  fora  =  1.0  and  7  -  5/3. 
Details  of  this  analysis  and  calculated  results  are  given  in  Subsection 
III-B.  This  work  represents  extension  of  our  previous  linearized  model  for  a 
single  component  system  to  include  nonlinear  effects.  In  Subsection  III-C, 
nonlinear  balance  equations  are  derived  considering  multiple  species  and  the 
bounding  cases  of  frozen  and  equilibrium  chemistry  in  the  Knudsen  layer.  This 
work  represents  extension  of  our  previous  linearized  multiple  species  model  to 
include  nonlinear  effects.  Calculated  results  for  this  model  also  are  given. 
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B.  SINGLE  COMPONENT  NONLINEAR  MODEL 


1.  Problem  Formulation 


Following  the  approach  of  Anisimov^,  as  extended  by  Knight*®  to 
include  rotational  degrees  of  freedom  of  the  molecules,  the  conservation  equa 
tions  of  mass  momentum  and  energy  for  the  Knudsen  layer  shown  in  Fig.  7  are 


pu  =  «ps 


merfc(m) 


(III-l) 


p(  S  +  RT) 


WTs 


(2-08  0PRT 


(m2  +  i )  erfc(m) 


(III-2) 
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RT  +  £ 


2RT$  +  G(y)R(T$ 


+O0pRT 


erfc(m)  -  (m2  +  2)e“m 
(III- 


where  G(y)  =  (5-3y)/(2(y-l)).  For  a  vaporization  coefficient  a  of  unity, 
these  equations  reduce  identically  to  those  given  by  Knight.  The  velocity 
distribution  functions  assumed  in  obtaining  these  equations  are: 


For  the  continuum  flow  at  the  edge  of  the  Knudsen  layer 


f 


3 


P(2jtRT)'3/2  exp 


2RT 


(II I -4a) 


For  the  saturated  vapor  (subscript  s)  molecules  ente',i.’io  the  vapor  phase  from 
the  solid  or  liquid 


-3/?  £ 2  ,  yj2  +  > 2 

fx  =  />s  (27rRTs)  j/^  exp  —  7RT  —  ~  '  j  •  ^  >0  (HI-4b) 

For  the  molecules  entering  the  solid  phase  or  reflecting  back  into  the  vapor 
phase  from  the  vapor  phase 

f2  s0f3  »  £  <  0  (II I -4c ) 

The  contributions  of  the  velocity  distribution  functions  fp  f^,  and  f^ 
to  corresponding  fluxes  J^,  J2,  and  in  the  conservation  equations  at 
the  interface  are  illustrated  in  Fig.  8  for  arbitrary  values  of  the  vaporiza¬ 
tion  and  condensation  coefficients.  In  obtaining  Eqs.  (III-l)  through 
(III -3) ,  it  has  been  assumed  that  the  condensation  coefficient  o»c  is  identi¬ 
cally  equal  to  the  vaporization  coefficient ay,  since  this  is  known  to  hold 
near  equilibrium;  thus,  c*  =  Ct  =  . 

w  V 

Equations  (III-l)  through  ( 1 1 1 -3 )  may  be  algebraically  manipulated  to 
obtain  analytical  expressions  for  T/Ts,  p/ps,  and  f}  as  functions  of  m  = 
u//  2RT  or  Mach  number  M  *  J (2/y) m.  Carrying  out  the  somewhat  tedious 
algebra,  we  obtain 
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Vi  j2  (I*“cu2 


Fig.  8.  Fluxes  J2,  and  J3  in  the  Knudsen  Layer 
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where 

<p(ce,m)  =  1  -  [m2  -  i r*^2  m(m2  +  ^  erf(m)em  j 

A(y,a,m)  =  5  +  4G(y)  -  r(a,m)  +  2m2  £l  -  ^(a,m)j 
B(y,or;m)  =  nl/2  ^  m?(of,m) 

It  can  be  easily  shown  that  for  a  =  1  and  y  =  5/3,  Eqs.  (III-5)  reduce  to 
those  given  by  Knight.  Note  that  although  Knight*®  wrote  his  conservation 
equations  for  arbitrary  values  of  y,  his  expressions  for  T/Ts,  p/ps>  and  fi 
are  valid  only  for  y*  5/3. 
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2. 


Calculated  Results 


For  specified  values  of  the  vaporization  coefficient  a  and  the  ratio  of 
specific  heats  7,  Eqs.  ( 1 1 1 -5 )  may  be  readily  solved  to  obtain  T/Ts  and  p/p$  as 
functions  of  the  Mach  number  at  the  edge  of  the  Knudsen  layer  M  =  (  \fTfy) m.  The 
pressure  ratio  p/p$  and  the  mass  flux  ratio  pu/  [  p$  J  RTs/27r ]  then  may  be 
obtained  from 


E_ 

Ps 


( 1 1 1 -6 a ) 


( I II-6b) 


To  a  first  approximation,  the  carbon  vapor  phase  may  be  assumed  to  con¬ 
sist  of  species  C-j  molecules  only.  Then,  for  ideal  gas  behavior  the  7  is 
equal  to  9/7.  Calculated  temperature,  density,  pressure,  and  mass  flux  ratios 
as  a  function  of  Mach  number  for  7  =  9/7  and  vaporization  coefficients  of  1.0, 
0.6,  and  0.2  are  shown  in  Fig.  9.  The  effect  of  vaporization  coefficient  on 
the  temperature  ratio  is  minimal.  However,  much  larger  effects  are  seen  in 
the  density,  pressure,  and  mass  flux  ratios.  For  these  latter  ratios,  a 
decrease  in  the  vaporization  coefficient  from  1.0  to  0.2  significantly 
decreases  the  Knudsen  layer  edge  quantity  relative  to  the  reference  saturation 
value  of  that  quantity.  Particularly  noticeable  are  the  large  drops  in  the 
density  and  pressure  ratios  for  a  =  0.2  between  Knudsen  layer  edge  Mach  num¬ 
bers  M  of  0  and  0.1  and  the  very  flat  mass  flux  ratio,  again  for  a  -  0.2, 
above  M  =  0.4. 


The  effect  of  the  ratio  of  specific  heats  7  for  unit  vaporization  coef¬ 
ficient  a  are  illustrated  in  Fig.  10.  The  effect  is  relatively  weak  for  all 
the  ratios,  especially  the  density  ratio.  In  Fig.  11,  predicted  results  from 
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Fig.  9.  Effect  of  Vaporization  Coefficient  a  on  Jump  Conditions 
Across  Knudsen  Layer 
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Fig.  10.  Effect  of  7  on  Jump  Conditions  Across 
Knudsen  Layer 
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near  Models 
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our  earlier  linearized  model^  and  the  present  nonlinear  model  for  y  -  9/7 
and  a  =  1.0  are  compared.  The  temperature  ratio  calculated  using  the  linear 
model  lies  below  that  predicted  by  the  nonlinear  model,  but  the  error  is  less 
than  eight  percent  even  for  a  Knudsen  layer  edge  Mach  number  M  of  0.5.  The 
pressure  ratio  predicted  by  the  linear  model  also  lies  below  that  of  the  non¬ 
linear  model.  However,  in  this  case  the  error  is  about  16  percent  for  M  =  0.3 
and  increases  dramatically  as  M  increases  further. 

3.  Application  to  Carbon  Nonequilibrium  Phase  Change 

The  results  presented  above  may  be  easily  applied  to  carbon  nonequili¬ 
brium  phase  change  analyses  if  it  is  assumed  that  the  carbon  vapor  is  made  up 
of  a  single  component,  e.g.,  carbon  species  C^.  Alternatively,  one  could 
assign  gross  properties  to  the  vapor,  i.e.,  average  molecular  weight,  mean 
values  of  heat  of  formation,  free  energy,  vaporization  coefficient,  and  so 
forth. 


The  calculations  are  closed  in  the  same  way  previously  discussed  using 
the  linearized  model^’^.  Equations  (lll-5a)  through  (III-6b)  or  Figs.  9  and 
10  provide  the  jump  conditions  across  the  Knudsen  layer,  i.e.,  the  ratio  of 
Knudsen  layer  edge  quantities  to  those  in  the  saturated  vapor  at  the  surface 
temperature,  as  a  function  of  Mach  number. 

To  calculate  the  mass  loss  rate  m  =pu  for  specified  values  of  the 
ambient  pressure  pg  and  input  radiation  heat  flux  IQ,  the  above  equations 
must  be  solved  together  with  the  surface  energy  balance  equation.  For 
steady-state  ablation,  this  equation  is 


I0  -  m  [h(T,p)  -  CpTb]  -<7€Ts4  =0  ( 1 1 1-7 ) 
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Assuming  that  the  carbon  vapor  is  made  up  of  species  molecules, 

the  mass  loss  rate  m  can  be  expressed  using  Eq.  ( 1 1 1-6b)  solely  as  a  function 

of  surface  temperature  T$  and  Mach  number,  i.e., 

ps  <y 

»  *  266  jjj-  (M)  ---I/!  ( II 1-8) 

2 

For  m  in  g/cm  -sec,  the  carbon  species  vapor  pressure  p$  should  be 
in  atmospheres  and  the  surface  temperature  T$  in  degrees  Kelvin. 

Calculation  of  m  thus  reduces  to  the  following  iterative  procedure: 

a.  Guess  Ts. 

b.  Calculate  ps(Ts)  and  p/ps  =  pe/p$. 

c.  Determine  the  Mach  number  from  Fig.  9  or  from  Eqs.  (II I-6a ) , 

(III-5a),  and  ( II I-5b) . 

d.  For  guessed  Ts  and  calculated  M,  determine  m  from  Eq.  (Ill -8 ) - 

e.  For  temperature  T  obtained  from  Eq.  (III-5a),  calculate  the 
enthalpy  of  the  vapor  at  the  edge  of  the  Knudsen  layer  h  = 
h(T,pe).  Using  this  value  of  h,  check  the  surface  energy 
balance,  Eq.  (III-7),  for  the  specified  value  of  ’0. 

f.  If  Eq.  ( 1 1 1 -7 )  is  not  satisfied,  guess  a  new  value  of  Ts  and 
return  to  step  a. 

Calculated  results  from  Eq.  (Ill -8 )  for  Mach  numbers  of  0.2,  0.4,  0.6, 
0.8,  and  1.0  are  shown  in  Fig.  12.  Thus,  utilizing  Figs.  9  through  12,  the 
iterative  solution  procedure  above  can  easily  be  carried  out.  Alternatively, 
the  equations  discussed  above  can  be  used  to  calculate  lines  of  constant  I 

0 

and  constant  pg  in  the  m  versus  1/T$  plane.  Then,  the  entire  iterative 
procedure  discussed  above  is  replaced  by  a  single  figure.  The  results  of  such 
calculations  for  the  present  model  are  shown  in  Fig.  13.  From  this  figure  for 
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specified  IQ  and  pg,  the  mass  loss  rate  m  and  the  surface  temperature  T$ 
can  be  obtained  directly  by  interpolation  between  plotted  values  of  I  and 

Pe‘ 

C.  MULTIPLE  COMPONENT  NONLINEAR  MODEL 

1.  Problem  Formulation  -  Frozen  Chemistry 

The  mass,  momentum,  and  energy  conservations  equations  ( 1 1 1 -1 )  through 
(III-3),  when  generalized  to  allow  for  more  than  one  species  in  the  vapor 
phase,  become  for  n  species 
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where  G(yi)  =  (5  -  ly^)/  2(^-1)  ,  subscript  i  designates  the  species  i 
values  of  the  quantity,  and  the  F  functions  are  given  by 


F*  =  7r1^mierfc(mi)  -  e  i 


(III-12a) 


Fi  s  (mi2  +  \  )  erfc(mi)  ‘  'ni7r"1/2  e 


-1/2  ~mi 


( 1 1 1 -12b) 


Fi  =  mi  ^mi  +  |^1/2erfc(mi) 


/2RT 


2  \  'mi  1 

m‘  +  2  j  e  1  (III-12C) 


( 1 1 1 - 12d ) 


In  order  to  solve  these  equations,  additional  information  is  required. 

A  portion  of  this  information  is  obtained  from  the  chemistry  of  the  system. 
Here  we  assume  frozen  chemistry,  i.e.,  no  chemical  reactions  at  all  take  place 
between  the  different  carbon  species.  Since  no  reactions  occur,  the  ratio  of 
the  species  i  mass  fraction  to  the  total  mass  fraction  must  be  the  same 
as  the  ratio  of  the  mass  flux  of  species  i  leaving  the  surface  to  the  total 
mass  flux.  In  other  words, 


Ki  _mi 

r "  r 


(III-13a) 


Substituting  K*  *  p^/p,  =  p^  and  m  =  pu,  Eq.  (III-13a)  becomes 


uf  s  u  (i*l,2,...n) 


(III-13b) 
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Additional  equations  required  to  close  the  system  are 


n 

p  -  £  pi 

i*l 


(II I -14a) 
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(II I -14d) 


and  the  surface  energy  balance  equation  (III-7),  where  m  and  h  are  now  given 
by 
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"■ s  pu  ■  2  Piui 
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(III-14e) 
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For  specified  IQ  and  pg,  Eqs.  ( 1 1 1 -9 )  through  (II 1-14 )  may  be 
solved  simultaneously  using  an  iterative  procedure  and  guessing  the  tempera¬ 
ture  Ts  as  described  in  Subsection  III-B.  For  the  guessed  temperature,  the 
carbon  species  vapor  pressure  p$  for  i  =  1,2. ..n  may  be  calculated  from 
thermochemical  data  and  subsequently  used  to  determine  the  . 

2.  Problem  Formulation  -  Equilibrium  Case 

Equilibrium  chemistry  in  the  Knudsen  layer  represents  the  opposite 
bounding  case  to  that  of  frozen  chemistry,  i.e.,  chemical  reactions  take  place 
between  all  the  carbon  species  at  rates  sufficient  to  ensure  local  chemical 
equilibrium.  Thus,  in  the  system  of  equations  given  by  Eqs.  ( 1 1 1 -9 )  through 
(III-14),  it  is  necessary  to  replace  those  equations  limited  to  frozen 
chemistry  by  alternate  equations  applicable  to  equilibrium  chemistry. 

When  chemical  reactions  take  place,  the  n  equations  given  by  Eq. 
(III-13b)  no  longer  are  valid.  However,  if  chemical  equilibrium  is  achieved 
at  the  edge  of  the  Knudsen  layer,  the  species  partial  pressures  p.  are  given 

by 


prm 


P-  (*0  i  -  2,3, .  • .  n 
i 


(III -15  a } 


where  T  is  the  temperature  at  the  edge  of  the  Knudsen  layer.  The  nth  inde¬ 
pendent  equation  is 
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A  basic  objective  of  this  work  is  to  obtain  carbon  species  thermo¬ 
chemical  data  information  from  laser  heating  experiments.  For  specified 
thermochemical  data  which  define  the  carbon  species  vapor  pressures  p$ 
and  assumed  vaporization  coefficients  a..,  the  equations  described  abov^ 
allow  the  mass  loss  rate  m  and  the  surface  temperature  T$  to  be  calculated 
as  a  function  of  the  ambient  pressure  pg  and  the  laser  radiation  heat  flux 
IQ.  The  results  of  calculations  using  a  somewhat  simplified  set  of  Eqs. 
(III-9)  through  (III-ll),  assuming  frozen  chemistry  in  the  Knudsen  layer, 
JANNAF  thermochemical  data,  and  unit  vaporization  coefficients,  are  shown  in 
Fig.  14. 

Comparison  of  these  predicted  results  with  laser-heated  carbon  ablation 

data  should  allow  one  to  determine  the  validity  of  the  JANNAF  data  and  the 

unit  vaporization  coefficient  assumption.  However,  as  discussed  in  our 

4  5 

earlier  work  with  linearized  models  ’  ,  a  uniqueness  problem  arises;  i.e., 
the  measured  mass  loss  rate  and  surface  temperature  can  be  predicted  using 
JANNAF  thermochemical  data  (which  give  relatively  low  carbon  total  vapor  pres¬ 
sure)  with  unit  vaporization  coefficients,  as  in  Fig.  14,  or  by  using  other 
thermochemical  data^’^,  which  give  carbon  total  vapor  pressures  up  to  an 
order  of  magnitude  larger  than  the  JANNAF  data,  with  reduced  values  of  the 
vaporization  coefficients. 


15  - 

Kratsch,  K.M.,  et  al.,  "Graphite  Ablation  in  High  Pressure  Environments," 
AIAA  Paper  68-1153,  1968. 

16 

Dolton,  T.A.,  H.E.  Goldstein,  and  R.E.  Mauer,  "Thermodynamic  Performance  of 
Carbon  in  Hyperthermal  Environments,"  Progress  in  Astronautics  and  Aero- 
nautics:  Thermal  Design  Principles  of  Spacecraft  and  Entry  Bodies,  21,  J.T. 
Bevans  (ed.),'  ATAAT  NewTorlc,  1959,  pp.  159-201'. - - 


In  an  attempt  to  resolve  the  uniqueness  problem,  sensitivity  calcula¬ 
tions  were  carried  out  to  determine  if  experimental  mass  loss  and  surface  tem¬ 
perature  data  could  be  taken  in  such  a  way  as  to  uniquely  determine  the  vapor 
pressure  (and  thus  the  thermochemical  data)  and  the  vaporization  coefficients 
from  the  data.  Mass  loss  rate  m  and  surface  temperature  Ts  were  treated  as 
dependent  variables  to  be  measured.  The  main  question  is,  How  should  the 
ambient  pressure  pg  and  the  laser  radiation  heat  flux  I  be  varied  in 
order  to  obtain  unique  information  from  the  measured  m  and  T$? 

An  early  conclusion,  in  line  with  previous  work  using  the  linearized 
nonequilibrium  phase  change  model,  was  that  if  all  the  data  are  taken  such 
that  the  outward  flow  at  the  surface  is  sonic  (solid  line  on  Fig.  14),  then 
uniqueness  cannot  be  established.  For  instance,  suppose  we  had  experimental 
data  represented  by  the  four  »'s  on  Fig.  14.  We  see  that  the  JANNAF  thermo¬ 
chemical  data  with  unit  vaporization  coefficients  predict  the  "data."  How¬ 
ever,  other  thermochemical  data  which  give  higher  total  vapor  pressure  also 
predict  the  "data"  using  smaller  vaporization  coefficients. 

To  eliminate  this  impasse,  sensitivity  calculations  were  made  for 
potential  experimental  data  in  the  subsonic  regime,  i.e.,  the  dashed  lines 
region  of  Fig.  14.  Two  chemistries  were  used:  JANNAF  and  Lee  and 
Sanborn^.  The  Lee  and  Sanborn  thermochemical  data  were  selected  to  compare 
with  JANNAF  because  they  are  more  recent  than  others  we  have  used  pre- 
viously^’^’^  and  give  carbon  total  vapor  pressures  only  2.5  to  3  times 
higher  than  JANNAF.  In  comparison,  thermochemical  data  of  Refs.  15  and  16 
result  in  predicted  carbon  total  vapor  pressures  an  order  of  magnitude  above 
JANNAF  for  temperatures  greater  than  4000°K.  Thus,  if  a  method  of  obtaining 
experimental  data  could  be  identified  which  would  allow  one  to  choose  uniquely 
between  JANNAF  and  Lee  and  Sanborn  thermochemical  property  data,  then  the  same 
method  surely  could  be  used  to  consider  other  thermochemical  data  which  differ 
from  JANNAF  to  a  much  greater  degree  than  that  of  Lee  and  Sanborn. 

17 - 

Lee,  E.L.  and  R.H.  Sanborn,  "Extended  and  Improved  Thermal  Functions  for  the 
Gaseous  Carbon  Species  C1-C7  from  298  to  1000  K,"  High  Temperature 
Science,  5  1973,  pp.  438-453. 
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It  was  found  that  for  an  assumed  carbon  species  C,  vaporization 

i  1  1  O 

coefficient  of  0.73,  and  also  assuming  a.  =  (a^)  ,  i  =  2,3,...n  , 

the  predicted  results  using  Lee  and  Sanborn  chemistry  matched  those  calculated 
using  JANNAF  chemistry  with  =  1.0,  and  both  calculations  matched  the 
sonic  line  ''data"  represented  by  the  four  *'s  on  Fig.  14. 

It  was  assumed  for  the  calculated  results  shown  in  Fig.  15  that  experi¬ 
mental  data  could  be  taken  in  such  a  manner  that  I  and  p„  varied  so  as  to 

o  e 

keep  the  surface  temperature  constant.  The  four  *'s  on  Fig.  15  correspond  to 
those  on  Fig.  14.  For  each  of  the  four  constant  values  of  surface  tempera¬ 
ture,  the  external  pressure  pg  was  increased  and  the  laser  radiation  heat 

flux  I  was  decreased  so  as  to  keep  T  constant.  The  calculated  results 
0  s 

show  very  significant  differences  in  behavior  between  the  JANNAF  chemistry  and 

Lee  and  Sanborn  chemistry  predictions.  Even  though  both  chemistries  match  the 

hypothetical  sonic  flow  data  represented  by  the  ^'s,  the  external  pressure 

increase  to  cause  an  order  of  magnitude  decrease  in  m  is  about  a  factor  of 

★ 

four  above  the  sonic  pressure  value  pg,  if  the  JANNAF  thermochemical  data 
are  correct,  whereas  a  factor  of  eight-nine  increase  is  required  if  the  Lee 
and  Sanborn  chemistry  is  correct.  This  Is  excellent  sensitivity,  but  data  are 
required  at  constant  temperature  which  is  difficult  to  obtain. 

It  is  much  easier  to  obtain  data  at  constant  ambient  pressure.  Similar 
calculations  for  these  two  chemistries  varying  the  laser  heat  flux  I  but 
keeping  ambient  pressure  constant  are  shown  in  Fig.  16.  Along  the  sonic  line, 
both  chemistries  match  the  four  "data"  points  with  different  vaporization 
coefficients  as  before.  However,  if  additional  data  could  be  obtained  at  con¬ 
stant  pressure,  it  should  be  possible  to  select  JANUAF  or  Lee  and  Sanborn 
chemistry,  especially  for  pg  =  1  atm  where  the  predicted  surface  tempera¬ 
tures  T$  for  the  nearly  vertical  portions  of  the  curve  differ  by  140°K  for 
the  two  chemistries. 

is - 

Rosenblatt,  G.,  Personal  Conmunication,  Department  of  Chemistry,  Pennsyl¬ 
vania  State  University,  University  Park,  Pa.,  Dec.  1979. 
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The  above  discussion  leads  to  a  conclusion  which  is  very  important  but 
has  not  been  previously  recognized.  It  is  illustrated  in  Fig.  17.  Whereas  in 
Fig.  16  two  chemistries  were  considered  with  different  vaporization  coef¬ 
ficients  for  each,  in  Fig.  17  JANNAF  frozen  chemistry  predicted  results  are 
shown  for  a  range  of  vaporization  coefficient  values.  It  is  obvious  from 
this  figure  that  the  nearly  vertical  rapid  rise  in  m  with  Ts  is  independent 
of  vaporization  coefficient.  Furthermore,  it  depends  only  on  the  vapor  pres¬ 
sure  and  the  ambient  pressure  which,  for  the  vertical  portion  of  the  curve, 
are  to  a  very  good  approximation  equal  to  one  another.  Thus,  we  conclude  from 
Figs.  16  and  17  that  a  rapid  rise  in  the  m  versus  1/T  curve,  from  experi¬ 
mental  data  at  constant  ambient  pressure  pg,  can  be  used  to  uniquely  deter¬ 
mine  vapor  pressure  at  the  experimentally  measured  value  of  Ts«  The  same 
data  as  well  as  other  data  along  the  sonic  line  could  then  be  used  to  deter¬ 
mine  vaporization  coefficients  from  measured  mass  loss  rate,  m,  values. 


IV.  CONTINUOUS  WAVE  LASER  DATA 


A.  MATERIALS.  EXPERIMENTAL  FACILITY.  TEST  APPARATUS. 

AND  INSTRUMENTATION 

1.  Materials 

Two  types  of  carbon  were  used  in  both  the  continuous  wave  (CW)  laser 
tests  reported  here  and  the  pulsed  laser  tests  discussed  in  Section  V.  For 
both  carbon  types,  samples  from  the  same  lot  numbers  were  used  in  the  CW  and 
pulsed  testing  to  ensure  meaningful  data  comparisons.  The  two  carbon  types 
tested  were  pyrolytic  graphite  (PG)  (Super  Temp  Company,  continuously 
nucleated)  and  Graphnol  (produced  by  Great  Lakes  Carbon  Research  Corporation 
for  the  Naval  Surface  Weapons  Center).  Cylindrical  samples  of  both  types  of 
carbon  with  dimensions  as  shown  in  Fig.  18  were  fabricated.  The  hole  in  one 
end  of  the  sample  facilitated  spinning  the  sample  by  mounting  it  on  a  mandrel 
which  could  be  rotated  about  its  axis  at  up  to  40,000  rpm.  For  the  PG,  the 
cylindrical  surface  was  a-face  and  the  flat  ends  were  c-face.  Slotted 
Graphnol  samples  were  prepared  when  it  became  apparent  from  initial  test  data 
on  the  nonslotted  Graphnol  samples  that  significant  convective  and  radiative 
heat  loss  was  occurring  due  to  conduction  toward  the  mandrel -support  end  of 
the  sample. 

2.  Experimental  Facility 

The  experiments  were  carried  out  using  the  CO2  electric  discharge 
laser  located  at  the  Philco-Ford  facility  at  Newport  Beach,  California.  This 
laser  had  been  used  previously  on  four  other  occasions  to  carry  out  similar 
testing.  A  sumnary  of  the  results  of  these  tests  recently  has  been  prepared 

1  q 

by  Whittaker,  et  al.  .  The  maximum  output  power  in  the  current  tests  was 

13 - 

Whittaker,  A.G.,  et  al.,  "Carbon  Vapor  Pressure  in  the  Range  3450-4500  K  and 
Evidence  for  Melting  at  ~3800  K,"  to  be  submitted  to  High  Temperature  Science. 
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6.4  kW.  The  10.6  /im  beam  was  approximately  8  cm  in  diameter  with  a  flat  (top 

hat)  power  distribution  across  the  beam.  Run  times  varied  from  2  to  12  sec, 

depending  upon  the  power  density  on  target  as  related  to  sample  survival  time 

and  the  specific  objectives  of  the  run.  Beam  power  was  measured  by  reflecting 

part  of  the  beam  to  a  calorimeter.  For  a  given  laser  power,  power  density  on 

the  target  could  be  varied  over  a  significant  range  by  using  a  variable  aper- 
20 

ture  centered  in  the  beam  in  front  of  the  calorimeter.  Estimated  time  to 
full  laser  power  was  1  to  1.5  sec  from  test  initiation.  To  compensate  for 
this,  a  quick -acting  reflective  shutter  was  available  and  utilized  in  a  por¬ 
tion  of  the  tests.  By  delaying  the  shutter  to  open  nominally  2  sec  after  test 
initiation,  full  laser  power  on  the  sample  from  the  initial  exposure  could  be 
assured. 

3.  Test  Apparatus  and  Instrumentation 

A  schematic  diagram  of  the  test  apparatus  is  shown  in  Fig.  19.  The 
pressurizable  test  chamber  was  made  of  brass  tubing  about  50  cm  long  with  an 
inside  diameter  of  just  over  50  cm.  The  end  plate,  viewing  window,  and 
smaller  NaCl  laser  windows  were  sealed  with  0-ring  seals.  The  focusing  lens 
reduced  the  beam  diameter  from  the  nominal  8  cm  output  diameter  received  from 
the  laser  to  a  3.5  iron  spot  diameter  on  the  rotating  cylindrical  carbon  sample. 

The  test  chamber  could  be  pressurized  from  0.1  to  800  torr  with  Ar, 

Oj,  or  arbitrary  mixtures  of  these  gases.  For  tests  conducted  at  one  torr, 
the  chamber  was  operated  in  an  open  configuration  with  a  continuous  flow  of 
the  pressurizing  gas  into  and  out  of  the  chamber.  For  higher  chamber  pres¬ 
sures,  a  closed  configuration  was  used.  Pressures  up  to  10.0  torr  were  mea¬ 
sured  using  a  Baratron  gage,  while  pressures  in  the  10  to  800  torr  range  were 
determined  using  a  25.4  cm  Heise  gage. 

to - 

Whittaker,  A.G.,  et  al.,  "System  Employing  Laser  Heating  for  the  Measurement 
of  High-temperature  Properties  of  Materials  over  a  Wide  Pressure  Range,"  Re- 
view  of  Scientific  Instruments,  48(6),  June  1972,  pp.  632-639. 
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Surface  temperatures  of  the  carbon  samples  were  measured  by  an  Ircon 
optical  pyrometer  which  viewed  the  rotating  sample  30  deg  away  from  the  laser 
beam  axis  in  the  direction  of  sample  rotation.  The  area  viewed  by  the  pyro¬ 
meter  was  a  circular  spot  approximately  1  nm  in  diameter.  Spectroscopic  data 
were  obtained  from  the  carbon  gas  plume  by  means  of  the  quartz  viewing  window 
at  the  top  of  the  chamber.  The  spectrograph  was  focused  at  a  point  2.5  mm 
from  the  sample  surface  on  the  axis  of  the  laser  beam.  For  selected  runs,  a 
high-speed  (Hycam)  camera  was  used  to  obtain  motion  pictures  (7000  frames  per 
second)  through  the  front  viewing  window,  and  particles  from  the  rotating 
sample  were  caught  in  a  tray  placed  in  the  chamber. 

An  eight-channel  Visicorder  was  used  to  record  data  during  each  test 
run.  Signals  from  the  laser  power  calorimeter,  pyrometer,  Baratron  pressure 
gage,  spectrograph  shutter,  and  the  time-delay  reflective  shutter  were  re¬ 
corded  in  time  sequence.  A  more  detailed  description  of  a  similar,  but  smal¬ 
ler,  test  chamber  with  the  same  instrumentation  is  given  by  Whittaker,  et 


B.  TEST  PLANNING 


Final  planning  for  this  test  series  was  carried  out  after  the  multiple 
component,  nonlinear  Knudsen  layer  model  described  in  Section  III  had  been 
developed  utilizing,  in  particular,  results  of  the  sensitivity  analyses.  As 
discussed  in  Subsection  III-C,  experimental  data  obtained  under  conditions 
(laser  power  and  ambient  pressure)  for  which  the  flow  in  the  laser-induced 
plume  away  from  the  ablating  surface  is  subsonic  are  particularly  useful. 

This  is  because  appropriate  data  of  this  type  can  be  used  to  uniquely  deter¬ 
mine  vapor  pressure  and  vaporization  coefficients. 

Consequently,  the  test  matrix  for  these  experiments  was  planned  using 
the  following  guidelines: 

1.  Obtain  almost  all  data  in  the  subsonic  plume  region. 

2.  Locate  experimentally  the  nearly  vertical  portions  of  the  m 
versus  1/TS  (see  Fig.  14)  curves  by  conducting  tests  at 
constant  ambient  pressure  for  a  series  of  laser  power  levels. 
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3. 


4. 


Locate  the  limiting  sonic  line  in  the  m  versus  1/TS  plane  by 

a.  Conducting  tests  at  several  moderate  to  high -laser 
power  levels  for  an  ambient  pressure  of  one  torr  (thus 
assuring  sonic  flow  at  the  Knudsen  layer  edge). 

b.  Conducting  tests  at  constant  laser  power  and  reducing 
the  chamber  pressure  incrementally  until  doing  so  no 
longer  reduces  the  surface  temperature,  i.e.,  the  flow 
has  choked  and  is  therefore  sonic  at  the  Knudsen  layer 
edge. 

Obtain  mass  loss  rate  m  for  all  tests  by  weighing  before  and 

after  the  tests. 


5.  Select  run  conditions  for  constant  ambient  pressure  tests  (2, 
above)  and  constant  laser  power  tests  (3,  above)  in  such  a  way 
as  to  obtain,  as  nearly  as  possible,  experimental  data  at 
constant  surface  temperature. 


C.  DATA  SUMMARY 


In  line  with  the  above  guidelines,  the  test  matrices  as  carried  out  are 
sumnarized  in  data  Tables  1  and  2  for  Graphnol  and  PG,  respectively.  Since 
laser  power  and  ambient  chamber  pressure  are  the  independent  variables  in 
these  experiments,  individual  values  of  these  parameters  for  each  data  run  are 
given  in  the  column  at  the  left  and  the  row  at  the  top,  respectively,  of  each 
table.  Within  the  box  corresponding  to  the  intersection  of  a  given  laser 
power  and  ambient  pressure,  the  sample  number,  length  of  the  run,  measured 
surface  temperature,  and  mass  loss  rate  are  given. 

Potentially  useful  data  were  obtained  on  a  total  of  about  85  samples  at 
laser  power  levels  ranging  from  1.0  to  6.4  kW  and  at  ambient  chamber  pressures 
of  0.5  to  800  torr.  Except  as  noted  in  the  tables,  for  all  pressure  levels 
above  one  torr  the  chamber  gas  was  an  equal  molar  mixture  of  Ar  and  Oj.  All 
tests  at  one  torr  or  below  were  pure  O2.  As  discussed  in  Ref.  19,  the  oxy¬ 
gen  was  necessary  in  order  to  eliminate  voluminous  soot  formation  which  can 
completely  obscure  the  pyrometer  data  and  significantly  reduce  the  power 
transmitted  to  the  sample  by  coating  the  laser  window  and  absorbing  laser 
energy  within  the  carbon  vapor  plume. 
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Table  1.  Graphnol  Data  (Continued) 


L3 


pe,  torr 

i  1 

_ 

38 

76 

152 

380 

_ 1 

Laser  Power 

! 

2.90 

81,  8  sec 
35650k 

4.64  mg/sec 

2.98  kW 

i 

j 

90C,  7  sec 

368  7°* 

17.1  mg/sec 

3.05  kW 

i 

i 

i 

54,  8  sec 
3648°* 

15.4  mg/sec 

5X  02 

2.98  kU 

95%  A r 

i 

84,  8  sec 

348  3°* 

1.71  mg/sec 

10%  02 

2.91  kU 

90%  Ar 

1 

85,  8  sec 
3516°* 

2.49  mg/sec 

23%  02 

3.41  kW 

77%  Ar 

47,  8  sec 

36  540k 

6.32  mg/sec 

3.19  kW 

1 

i 

89C,  6  sec 
3681°* 

16.9  mg/sec 

3.34  ktf 

82,  6  sec 
3609°* 

14.9  mg/sec 

1 

1  i 

3.41  kW 

57,  7  sec 

363  1  0k 

15.2  mg/sec 

3.63  kW 

48C,  5  sec 
37540K 

18.4  mg/sec 

3.99  kW 

49C,  5  sec 
37770* 

23.2  mg/sec 

4.07  kW 

22,  5  sec 
36000* 

4.21  kW 

69C,  5  sec 
37940* 

29.6  mg/sec 

4.29  kW 

8,  7  sec 
35870* 

18.2  mg/sec 

4.36  kW 

1,  8  sec 

368  7  0* 

17.7  mg/sec 

Table  1.  Graphnol  Data  (Concluded) 
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Pyrolytic  Graphite  Data 


Table  2.  Pyrolytic  Graphite 


U3  mq/sec 


An  equal  molar  mixture  of  Ar  and  was  selected  because  this  per¬ 
centage  of  was  the  minimum  amount  required  for  ambient  pressures  near  76 
torr,  and  it  was  felt  that  maintaining  a  constant  mixture  ratio  for  all  total 
pressures  was  desirable.  At  ambient  pressures  significantly  below  76  torr, 
C>2  mole  fractions  greater  than  50  percent  are  required.  Thus,  data  at  one 
torr  were  obtained  using  100  percent  02* 

D.  PRELIMINARY  DATA  ASSESSMENT 


The  tests  were  conducted  in  the  tenth  month  of  the  contract  year.  The 
time  since  then  has  been  devoted  mainly  to  reducing  the  data  to  obtain  mass 
loss,  mass  loss  rate,  and  sample  surface  temperatures  for  each  run.  This 
information  is  summarized  in  Tables  1  and  2  for  the  Graphnol  and  PG  samples, 
respectively.  In  terms  of  the  guidelines  discussed  in  Subsection  IV-8,  the 
scope  of  the  data  obtained  is  considered  to  be  very  good.  In  the  following, 
the  post-test  appearance  of  selected  samples  of  both  types  of  carbon  are  first 
discussed  in  order  to  put  planned  detailed  data  analysis  in  the  proper  per¬ 
spective.  Then,  a  preliminary  data  analysis  is  presented,  comparing  predicted 
results  with  experimental  data  with  the  objective  of  obtaining  carbon  vapor 
pressure  and  vaporization  coefficient  information. 

1.  Post-Test  Appearance  -  Graphnol  Samples 

The  cylindrical  carbon  samples  with  pretest  geometries  as  shown  in  Fig. 
18  were  irradiated  with  the  3.5  mm  diameter  laser  beam  focused  toward  the  end 
of  the  sample,  where  the  carbon  material  being  tested  extended  all  the  way  to 
the  cylinder  axis.  This  was  done  in  order  to  maximize  the  amount  of  material 
available  to  be  ablated  away,  as  well  as  to  minimize  the  energy  lost  by  heat 
conduction  into  the  spinner  mandrel.  Most  samples  were  spun  at  a  rate  of 
20,000  rpm. 
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For  the  Graphnol  samples,  this  resulted  in  a  relatively  symmetric 

groove  of  various  depths,  depending  upon  the  laser  power.  An  example  of  such 

a  groove  for  a  low  pressure  (one  torr)  run  of  8-sec  duration  with  a  peak  laser 

power  of  4.35  kW  is  shown  in  Fig.  20.  The  surface  of  the  ablated  carbon  in 

the  groove  is  relatively  smooth  as  a  result  of  the  fine  grain  structure  of  the 

Graphnol  sample.  Outside  the  groove  on  the  cylinder  surface,  a  wide  band  of 

vapor-deposited  carbon  with  a  grey  surface  appearance  is  seen.  This  material 

21  22 

may  be  one  of  the  carbynes  recently  discussed  by  Whittaker  .  The  black 
material  visible  on  top  of  the  sample  is  soot. 

Further  examples  of  the  post-test  appearance  of  Graphnol  samples  are 
shown  in  Fig.  21.  Each  of  the  three  samples  was  exposed  to  a  nominal  laser 
power  of  3  kW.  The  ambient  pressures  left  to  right  were  380  torr,  38  torr, 
and  one  torr.  The  width  of  the  vapor -deposited  grey  band  decreases  signifi¬ 
cantly  at  the  higher  ambient  pressures.  This  is  because  the  lateral  expansion 
of  the  laser-induced  plume  of  carbon  vapor  is  greatly  reduced  as  the  ambient 
pressure  is  increased.  Of  greatest  significance  are  the  measured  mass  losses 
as  given  on  the  figure.  Since  the  nominal  laser  power  and,  therefore,  the 
power  density  on  the  samples  was  the  same  for  all  three  samples,  all  of  our 
previous  linearized  models  as  well  as  the  multiple  species  nonlinear  model  of 
Subsection  III-C  would  predict  the  mass  loss  rate  for  the  three  samples  to  be 
about  the  same.  The  very  large  increase  in  mass  loss  rate  is  due  to  the  pre¬ 
sence  of  O2  in  the  chamber  atmosphere.  Thus,  in  the  detailed  data  analysis, 
which  will  be  proposed  as  a  part  of  next  year's  contract,  mass  loss  rates  for 

Graphnol  sanples  will  need  to  be  corrected  for  mass  loss  due  to  oxidation  of 
3  23 

solid  carbon  ’  before  vaporization  coefficient  information  can  be  ob¬ 
tained.  Most  of  the  mass  loss  by  oxidation  is  believed  to  occur  from 

n - 

Whittaker,  A.G.,  Personal  Conmunication,  The  Aerospace  Corporation,  El 
Segundo,  Calif. 

22 

Whittaker,  A.G. ,  "Carbon:  A  New  View  of  Its  High-Temperature  Behavior," 
Science,  20,  19  May  1978,  pp.  763-764. 

23 

Doak,  R.B.  and  J.R.  Baron,  "Activation  Studies  in  Graphite  Oxidation,"  Pro- 

?ress  in  Astronautics  and  Aeronautics:  Rarefied  Gas  Dynamics.  51(11),  J.L. 
otter  (ed.),  AIAA,  New  York,  1977,  pp.  621-634. 
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Ambient  Pressure:  1  torr 
Laser  Power:  4.36  kW 
Duration:  12  sec  (Sample  1) 


20.  Graphnol  Sample,  Long  Duration,  Low 
Ambient  Pressure  Run 


Ambient  Pressure: 
Laser  Power: 
Duration: 

Mass  Loss  Rate: 


380  torr 
3.05  kW 

8  sec  (Sample  54) 
13.4  mg/sec 


38  torr 
3.05  kW 

6  sec  (Sample  83) 
9.15  mg/sec 


1  torr 
2.98  kW 

8  sec  (Sample  79) 
4.9  mg/sec 


Fig.  21.  Effect  of  Ambient  Pressure  at  Constant 
Laser  Power,  Graphnol  Samples 


the  cylindrical  surface  away  from  the  laser  groove.  For  Graphnol,  this  entire 
surface  becomes  glowing  hot  during  each  run. 

Other  factors  which  somewhat  complicate  the  accurate  determination  of  mass 
loss  rates  for  the  Graphnol  samples,  but  at  the  same  time  supply  potentially 
very  useful  new  information  concerning  the  high  temperature  behavior  of  car¬ 
bon,  are  shown  in  Figs.  22  and  23.  The  sample  shown  in  Fig.  22  was  exposed  to 
a  high  power  density  for  a  short  time  with  full  power  on  the  sample  from  its 
initial  exposure  through  use  of  the  time-delayed  reflective  shutter.  The 
ambient  pressure  was  relatively  low.  This  resulted  in  the  unique  formation  of 
the  vapor-deposited  globules  barely  visible  on  top  of  the  sample  in  the  left 
picture  and  shown  in  plan  view  on  the  right.  In  the  past,  similar  globules 
have  been  interpreted  at  times  as  evidence  of  liquid  carbon.  This  is  highly 
unlikely  in  this  case,  as  the  maximum  temperature  reached  by  the  sample  was 
only  3582°K.  It  is  speculated  that  the  globules  may  be  dendrites^.  On 
the  cylindrical  surface,  a  thin  vapor-deposited  band  of  the  previously 
described  grey  material  can  be  seen  along  with  a  much  wider  band  of  soot. 

The  five  samples  in  the  pictures  of  Fig.  22  show  the  post-test  appear¬ 
ance  of  slotted  Graphnol  samples.  From  Table  1,  it  is  apparent  that  the  slot 
did  somewhat  reduce  energy  loss  due  to  conduction  in  that,  for  identical  laser 
power  and  amfcient  pressure,  a  slotted  Graphnol  sample  is  likely  to  be  20  to 
30°K  higher  in  surface  temperature.  However,  the  slot  also  resulted  in 
vapor  deposition  on  the  relatively  cold  surface  formed  by  the  slot  and  the 
main  portion  of  the  cylinder.  In  the  upper  picture,  the  carbon  vapor- 
deposited  ring  is  shown  detached  from  the  surface  on  which  it  formed  due  to 
differential  expansion  on  cooldown.  The  sample  on  the  right  did  not  form  a 
ring,  possibly  because  of  the  relatively  low  ambient  pressure.  The  three 
samples  in  the  lower  picture  each  have  detached  vapor -deposited  rings  as  well 
as  other  possibly  vapor-deposited  material  on  the  outside  lip  of  the  ablative 
mass  loss  groove. 


Pressure:  38  torr 


Graphnol  Sample  with  Vapor-Deposited  Globules 


5.23  kW 

5  sec  (Sample  47) 


4.36  kW 

5  sec  (Sample  46) 


380  torr 
3.05  kW 

8  sec  (Sample  88) 


800  torr 

kw 

~3.5  sec  (Sample  93) 


I! 


Ill 


Ambient  Pressure:  380  torr 
Laser  Power:  4.5  kW 
Duration:  5  sec  (Sample  92) 


Fig.  23.  Slotted  Graphnol  Samples 


2.  Post-Test  Appearance  -  Pyrolytic  Graphite  Samples 

In  contrast  to  the  Graphnol  samples,  the  ablation  groove  for  the  PG 
samples  always  showed  an  irregular  profile.  Examples  of  this  are  shown  in 
Fig.  24.  The  two  samples  each  saw  about  the  same  nominal  laser  power,  the 
difference  between  the  two  being  a  fairly  substantial  difference  in  ambient 
pressure.  The  low  pressure  (152  torr)  sample  on  the  left  shows  an  irregular 
groove  and  a  silver-grey  deposit  on  the  cylinder  surface  outside  the  groove 
which  extends  down  into  the  groove.  The  high  pressure  (760  torr)  sample  on 
the  right  has  an  irregular  groove  and  a  distinctly  different  surface  appear¬ 
ance.  The  mass  loss  rates  measured  for  these  two  samples  are  within  10  to  15 
percent  of  each  other,  as  would  be  expected,  and  in  contrast  to  the  Graphnol 
samples  shown  in  Fig.  21  and  discussed  previously. 

Thus,  even  though  the  material  does  not  respond  uniformly,  with  gouges 
appearing  between  graphite  planes,  the  very  low  thermal  conductivity  perpen¬ 
dicular  to  the  c-faces,  i.e.,  in  the  direction  parallel  to  the  cylinder  axis, 
localizes  the  heating  to  a  fairly  narrow  strip  not  too  much  wider  than  the 
diameter  of  the  laser  beam.  This  prevents  significant  mass  loss  by  oxidation 
of  solid  carbon  on  the  cylindrical  surface,  as  occurs  for  the  Graphnol 
sanples.  This  view  is  substantiated  by  the  immediate  post-test  appearance  of 
the  samples  just  seconds  after  the  laser  has  shut  down.  The  PG  samples  at 
that  time  are  glowing  hot  only  in  a  band  centered  on  the  deepest  part  of  the 
laser -produced  groove  in  the  sample.  In  contrast  to  this,  the  Graphnol 
samples  at  that  time  are  glowing  hot  over  the  whole  cylindrical  surface. 

PG  samples  exposed  to  a  higher  laser  power  and  different  pressures  are 
shown  in  Fig.  25.  In  this  case,  the  sample  on  the  left  showed  significantly 
less  mass  loss  (at  the  same  laser  power)  due  to  voluminous  soot  formation 
which  blackened  the  laser  window  and  reduced  the  laser  power  density  on  the 
sample.  A  fairly  thick  and  wide  band  of  vapor-deposited  soot  can  be  seen  on 
the  cylindrical  surface  below  the  silver-grey  band.  The  groove  on  this  sample 
is  highly  irregular. 
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Pressure:  152  torr  760  torr 

wer:  4.14  kW  3.92  kW 


Graphite  Samples,  Effect  of 
at  Constant  Laser  Power 


To  obtain  carbon  vapor  pressure  and  vaporization  coefficient  infor¬ 
mation  from  the  experimental  data,  the  variation  of  the  measured  mass  loss 
rates  and  surface  temperatures  is  to  be  compared  with  theoretical  predictions 
as  discussed  in  Subsection  III-C.  A  preliminary  comparison  of  this  type  is 
shown  in  Fig.  26.  The  mass  loss  rates  plotted  in  this  figure  are  those  tabu¬ 
lated  in  Tables  1  and  2.  They  were  obtained  simply  by  dividing  the  total  mea¬ 
sured  mass  loss  by  the  total  sample  exposure  time.  Thus,  no  corrections  for 
transients  near  time  zero,  mass  loss  by  oxidation,  or  material  vaporized  and 
subsequently  redeposited  from  the  vapor  have  been  made.  However,  the  overall 
slope  of  the  data  in  general  follows  the  slope  of  the  reference  JANNAF  line 
drawn  for  comparison. 

Graphnol  data  are  given  by  the  circles  and  PG  data  by  the  triangles. 
Filled  symbols  designate  run  conditions  believed  to  correspond  to  sonic  flow 
at  the  Knudsen  layer  edge,  while  open  symbols  designate  run  conditions  pre¬ 
dicted  to  give  subsonic  flow.  The  PG  data  tend  to  lie  below  the  data  for  the 
Graphnol,  most  likely  due  to  significant  mass  loss  by  oxidation  in  the  latter 
case.  To  obtain  vapor  pressure  and  vaporization  coefficient  information  from 
these  data,  corrections  for  oxidative  mass  loss,  transient  heating,  and  vapor 
deposition,  as  well  as  correlation  with  the  macroscopic  energy  balance  to 
relate  the  data  points  to  the  laser  power,  will  be  required.  We  plan  to  pro¬ 
ceed  with  the  work  necessary  to  accomplish  these  requirements. 


V.  PULSED  LASER  DATA 


The  two  carbon  types  tested  in  these  experiments  were  the  same  as  those 
described  in  Subsection  IV-A,  pertaining  to  the  CW  laser  tests.  For  each  car¬ 
bon  type,  samples  from  the  same  lot  numbers  were  used  in  the  CW  and  pulsed 

laser  testing.  Whereas  the  power  densities  on  the  carbon  samples  were  10  to 
2 

50  kW/cm  for  the  CW  laser  10.6 //m  radiation,  the  power  densities  in  the 
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pulsed  tests  discussed  here  were  200  kW/cm  up  to  4  MW/cm  for  the  1.06 fim 
neodymium-glass  laser  radiation.  The  corresponding  maximum  exposure  times 
were  about  10  sec  and  1  msec,  respectively.  The  data  discussed  here  were  ob¬ 
tained  during  the  past  year  as  part  of  a  cooperative  effort  with  the  NASA  Ames 

24 

Research  Center  and  have  been  documented  recently  by  Covington  . 

A.  CARBON  VAPORIZATION  EXPERIMENTS 


1.  Background 

The  laser,  experimental  technique,  and  apparatus  used  in  these  experi- 
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ments  are  described  by  Covington,  et  al.  .  The  method  utilized  by 

Covington  and  coworkers  seeks  to  relate  an  experimentally  derived  stagnation 
pressure  to  the  vapor  pressure  of  refractory  materials.  The  stagnation  pres¬ 
sure  is  determined  by  measuring  the  geometry  of  the  free-jet  expansions  from 
laser-vaporized  surfaces  and  employing  a  semi -empirical  relationship  which 
relates  the  measured  geometry  to  the  stagnation  pressure  of  the  flow  and  the 
ambient  pressure.  In  Ref.  25,  it  was  assumed  that  the  vapor  pressure  of  the 
material  is  equal  to  the  stagnation  pressure  derived  from  the  experimental 
data.  Carbon  vapor  pressures  determined  in  this  manner  were  shown  to  agree 
reasonably  well  with  those  obtained  from  JANNAF^  thermochemical  property 
data  at  measured  surface  temperatures. 

Covington,  A.,  "Ames  Free-jet  Laser  Vaporization  Experiments,"  Data  Report, 
NASA  Ames  Research  Center,  Mountain  View,  Calif.,  Sept.  1980. 
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Covington,  M.A.,  G.N.  Liu,  and  K.A.  Lincoln,  "Free-jet  Expansions  from 
Laser-vaporized  Planar  Surfaces,"  AIAA  Journal ,  15,  Aug.  1977,  pp.  1174-1179. 
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Additional  analysis  of  these  experiments  by  Baker  and  Crowell  indi¬ 
cated  possible  carbon  melting  at  about  4250°K,  and  also  showed  that  inter¬ 
pretation  of  the  data  in  terms  of  JANNAF  properties  is  not  unique,  i.e.,  other 
thermochemical  data^’^  also  are  able  to  predict  the  measured  data.  The 
primary  objectives  of  the  new  tests  carried  out  this  past  year  have  been  to 
uniquely  determine  the  correct  thermochemical  property  data  (vapor  pressure) 

and  to  establish  the  melt  temperature  and  carbon  species  vaporization  coeffi- 
2 

cients  .  Progress  toward  meeting  these  objectives  is  surtinarized  below. 

2.  New  Data  Summary 

New  pulsed  laser  data  obtained  this  year,  as  reported  in  Ref.  24,  are 
sunmarized  in  Figs.  27  through  29.  Over  100  new  data  points  were  obtained, 
with  approximately  two-thirds  of  these  being  for  PG  samples  and  the  remainder 
for  Graphnol.  A  few  new  data  points  also  were  obtained  for  ATJ-S  graphite 
samples. 

The  PG  data  are  shown  in  Fig.  27.  For  each  data  point,  the  peak  total 
(stagnation)  pressure,  as  deduced  from  the  measured  geometry  of  the  free-jet 
expansion,  is  plotted  versus  the  measured  peak  surface  temperature.  Experi¬ 
mental  data  obtained  for  Graphnol  samples  and  new  (1980)  and  old  (1976)  ATJ-S 
samples  are  given  in  Figs.  28  and  29.  In  Fig.  28,  the  total  pressures  derived 
from  the  free-jet  geometry  measurements  are  shown  as  a  function  of  peak  laser 
flux  (power  density).  'Reasonable  agreement  is  seen  between  the  new  and  the 
old  ATJ-S  graphite  data  points  as  well  as  between  all  three  different  types  of 

carbon  at  a  given  value  of  q_  .  The  divergence  of  the  PG  data  from  the 

m  dx  c  9 

ATJ-S  and  Graphnol  data  below  6  *  5  x  103  W/cnr  is  due  to  the  power 

m  aX 

density  being  too  low  to  bring  the  sample  up  to  the  steady-state  ablation  tem¬ 
perature  in  the  ATJ-S  and  Graphnol  cases,  whereas  the  much  lower  PG  thermal 
conductivity  (perpendicular  to  c-face)  resulted  in  the  PG  reaching  the  abla¬ 
tion  temperature  during  the  laser  pulse  time. 
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Fig.  29.  Peak  Surface  Temperature  versus  Peak 
Laser  Flux  for  Vaporization  in  Argon 


In  Fig.  29,  the  measured  surface  temperatures  for  the  data  points  in 

Fig.  28  are  shown  for  the  various  peak  laser  flux  levels.  The  large  d if - 

-1  4-1 

ferences  in  the  data  for  inverse  temperatures  (T  x  10  ,  K  )  greater 

than  2.3  again  are  due  to  thermal  conductivity  differences  which  result  in 

•  5  2 

different  measured  peak  temperatures  when  q  is  below  about  5  x  10  W/cm  . 

m  six 

However,  for  inverse  temperatures  less  than  2.3,  other  differences  not 
attributable  to  thermal  conductivity  are  seen.  First,  the  new  and  the  old 
ATJ-S  data  are  not  in  agreement.  Also,  the  new  ATJ-S  and  Graphnol  data  group 
together  but  are  not  in  agreement  with  the  PG  data.  A  preliminary  analysis  of 
these  results  is  given  in  Subsection  V-D. 

B.  Surface  Reflectance  Measurements  and  Spectroscopy 

Surface  reflectance  measurements,  using  an  auxiliary  laser,  have  been 
shown  to  give  a  good  indication  of  surface  melting  associated  with  pulsed 

pc 

laser  annealing  of  semiconductors  .  Significant  progress  has  been  made 

this  year  toward  using  this  technique  to  detect  the  presence  of  a  thin  melt 

layer  on  an  ablating  carbon  surface.  The  auxiliary  laser  and  associated 

detection  equipment  have  been  purchased  and  put  into  place,  and  initial  mea- 

27 

surements  are  expected  to  be  made  in  the  fall  of  1980  . 

A  detailed  evaluation  of  pulsed  laser  experimental  data  will  require 
carbon  species  concentration  measurements  in  addition  to  the  stagnation  pres¬ 
sure  and  surface  temperature  data  as  given  in  Figs.  27  through  29.  By  com¬ 
parison  of  experimentally  determined  carbon  species  ratios  with  predictions 
for  assumed  frozen  and  equilibrium  kinetics  in  the  Knudsen  layer  (see  Sub¬ 
section  III-C),  a  scaling  relationship  for  the  appropriate,  use  of  these 

25 

Auston,  D.H.,  et  al.,  "Time  Resolved  Reflectivity  of  Ion-implanted  Silicon 
During  Laser  Annealing,"  Appl .  Phys.  Lett.,  33,  1978,  pp.  437-440. 
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Covington,  M.A.,  Personal  Comnunication,  NASA  Ames  Research  Center,  Moffett 
Field,  Calif. 
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bounding  cases  for  the  prediction  of  carbon  vaporization  at  low  pressures 
under  high  thermal  radiation  fluxes  can  be  determined.  Toward  this  end,  ini¬ 
tial  spectroscopic  data  of  the  type  intended  to  be  used  to  determine  the 
desired  carbon  species  ratios  are  shown  in  Fig.  30  (from  Ref.  24).  The  mea¬ 
sured  relative  intensity  levels  of  the  carbon  species  C7  (Swing  band)  emis- 
sion  spectrum  (A  =  3850  -  4200  A)  and  the  C9  (Swan  band)  emission  spectrum 

o  L~ 

(\=  5050  -  5200  A)  are  compard  with  predicted  results  from  a  synthesized  band 

28 

spectra  computer  program  .  We  plan  to  obtain  additional  data  of  this  type 
in  the  coming  year. 

C.  Thulium  Vaporization  Experiments 

24 

These  experiments  were  undertaken  to  test  the  validity  of  the 

method  being  used  to  deduce  carbon  vapor  pressure  information  from  laser- 

induced  free-jet  expansions  for  a  material  which  (a)  vaporizes  as  a  monatomic 

gas,  (b)  is  very  likely  to  have  a  unit  vaporization  coefficient,  and  (c)  has 

established  thermochemical  data  (vapor  pressure  information).  Thus,  the  ori- 

14 

ginal  model  of  Anisimov  is  applicable  to  such  data  without  the  further 
complicating  factors  (see  Subsection  III-C)  of  internal  degrees  of  freedom  of 
the  molecules  (polyatomic  gases),  nonunit  vaporization  coefficients,  and 
several  polyatomic  species  in  the  vapor  with  potential  chemical  reactions, 
i.e. ,  the  carbon  system. 

The  element  thulium  was  selected  because  it  was  desirable  to  have  a 
material  with  a  high  melt  temperature  and  significant  vapor  pressure  at  tem¬ 
peratures  just  below  the  melting  point.  These  criteria  were  an  attempt  to 
ensure  visible  emission  from  the  laser-vaporized  plume  so  that  it  could  be 
photographed  to  determine  its  geometry,  sublimation  of  the  material  to  avoid 
the  complications  of  a  melt  layer  and,  finally,  sufficiently  high  mass  loss 
rates  to  have  a  measurable  mass  loss.  Thulium  melts  at  about  1800°K  and 

25 - - - 

Arnold,  J.E.,  et  al.,  "Line-by-line  Transport  Calculations  for  Jupiter  Entry 
Probes,"  Progress  in  Astronautics  and  Aeronautics:  Entry  Heating  and  Thermal 
Protection.  69,  W.  Olstad  (ed.j,  AIAA.  New  York,  1980,  pp.  52-82. 
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Fig.  30.  Comparison  of  Calculated  and  Experimental 
Spectra— C2  {  Av  =  0)  and  C3 


attains  a  vapor  pressure  greater  than  250  torr  before  it  melts.  It  was 

clearly  the  best  of  readily  available  materials  with  the  desired  properties. 

Calculations  and  later  experimental  data  revealed  that  with  the  pulsed 
laser  available,  the  thulium  could  not  be  brought  up  to  a  steady-state  abla¬ 
tion  temperature  in  the  pulse  time  available  at  a  power  density  (laser  flux) 
low  enough  to  not  melt  the  surface.  A  representative  picture  of  a  laser- 
induced  plume  from  the  surface  of  a  flat  thulium  test  sample  is  shown  in  Fig. 
31.  In  this  picture,  and  many  others  taken  at  different  laser  flux  levels, 

ambient  chamber  pressures,  and  camera  exposure  times,  the  Mach  disc  location 

needed  to  deduce  plume  total  pressures  could  not  be  seen.  Significant  absorp¬ 
tion  and  emission  from  induced  fluorescence  of  the  thulium  atoms  is  taking 
place,  but  this  obscures  any  details  of  the  free-jet  flow  patterns  and  also 
tends  to  overshadow  the  continuum  radiation  from  the  hot  vaporizing  surface 
which  is  used  to  determine  ablation  temperatures. 

While  the  above  results  did  not  give  the  desired  data  for  verifying  the 

24 

modeling  and  data  analysis  procedures,  additional  measurements  by  Lincoln 
with  a  time-of -flight  mass  spectrometer  verified  that  the  major  species  in  the 
laser-induced  plume  is  monatomic  thulium.  Thus,  thulium  still  represents  a 
good  candidate  material  for  model  verification  if  an  independent  method  of 
determining  the  mass  loss  rate  can  be  found.  One  possibility  would  be  to  con¬ 
duct  CW  laser  tests  as  described  in  Section  IV  and  to  weigh  the  samples  before 
and  after  exposure.  Approximate  calculations  indicate  that  lower  laser  beam 
power  densities  and  longer  run  times  would  facilitate  obtaining  data  below  the 
thulium  melt  temperature. 

D.  PRELIMINARY  DATA  ANALYSIS 

This  discussion  is  primarily  limited  to  the  new  PG  data  shown  in  Fig. 
27.  The  sensitivity  analysis  performed  with  the  nonlinear  multiple  carbon 
species  Knudsen  layer  model,  as  discussed  in  Subsection  III-C,  clearly  shows 


Peak  Laser  Flux: 
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Laser  Beam  Diameter: 
Beam  Incident  Angle: 
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Ambient  Pressure: 
Surface  Temperature: 


2.0  x  105  W/cm? 

~  700  // sec 

2.4  rrm 

250 

Argon 

3.1  torr 

>  2000°K 


Fig.  31.  Self-Luminous  Monatomic  Thulium  Vapor  Plume 
Created  by  Laser  Heating 
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that  carbon  vapor  pressure  and  vaporization  coefficient  information  cannot  be 
unambiguously  obtained  from  data  in  which  the  flow  at  the  edge  of  the  Knudsen 
layer  is  always  sonic.  All  the  data  in  Fig.  27  fit  into  this  category.  Thus, 
the  final  interpretation  of  these  data  will  need  to  be  done  in  conjunction 
with  the  Section  IV  CW  laser  data  which  do  not  have  this  limitation. 

It  is  important  and  significant,  however,  to  realize  that  the  data  in 
Fig.  27  cover  measured  surface  temperatures  from  about  3850°K  up  to  greater 
than  4700°K.  In  the  temperature  range  from  3850  to  4000°K,  overlapping 
data  from  both  the  CW  and  pulsed  laser  tests  are  available.  This  provides  the 
link  required  for  unique  determination  of  vapor  pressures  and  vaporization 
coefficients.  Since  most  experimentally  determined  melt  temperatures  lie  be¬ 
tween  4000  and  4400°K,  the  surface  reflectivity  measurements  to  be  performed 
this  fall  to  look  for  surface  melting  will  aid  considerably  in  the  detailed 
interpretation  of  data  up  to  4700°K.  It  is  possible  that  melting  phenomena 
are  associated  with  some  of  the  anomalies  of  the  Fig.  29  data. 

While  published  experimental  data  indicate  melting  at  least  at  tempera¬ 
tures  above  4400°K,  the  theoretical  calculations  of  Ref.  29  predict  a  melt 
temperature  of  4765°K.  If  we  assume  no  melting,  even  at  4700°K,  what  do 
the  Fig.  27  data  indicate?  First  of  all,  the  slope  of  the  data  is  not  paral¬ 
lel  either  to  the  JANNAF  or  the  Ref.  29  predicted  vapor  pressures.  From  the 
analysis  of  Subsection  III-C,  it  can  easily  be  shown  that  for  these  experi¬ 
ments  the  vapor  pressure  is  proportional  to  the  measured  stagnation  pressure. 
Thus,  the  data  shoulo  be  parallel  to  the  vapor  pressure  curves  unless  a 
strongly  temperature-dependent  vaporization  coefficient  alters  the  p°  versus 
1/T$  behavior.  An  example  of  such  behavior  predicted,  with  vaporization 
coefficient  a  strong  function  of  temperature,  is  shown  in  Fig.  32. 

79 - 

Leider,  H.R.,  O.H.  Krikorian,  and  D.A.  Young,  "Thermodynamic  Properties  of 
Carbon  Up  to  the  Critical  Point,"  Carbon,  11,  1973,  pp.  555-563. 
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Fig.  32.  Dependence  of  Mass  Loss  Rate  on 
Vaporization  Coefficient  and  Ts 


Even  with  the  carbon  species  vaporization  coefficient  decreased  to  its 
experimentally  determined  value  of  0.023  at  2700°K  and  increasing  to  1.0  at 
high  temperatures  (a  strong  temperature  dependence),  the  predicted  mass  loss 
rates  m  do  not  fall  off  at  lower  temperatures  as  rapidly  as  the  data  in 
Fig.  27  would  seem  to  indicate.  Thus,  other  possible  causes  of  this  behavior, 
including  chemical  reactions  or  condensation  occurring  in  the  plume  and  af¬ 
fecting  the  measured  Mach  disc  location  measurements,  are  being  investigated. 


VI.  CARBON  MELT  TEMPERATURE 


Controversy  concerning  the  melt  temperature  of  carbon  has  had  a  long 
30 

history  .  While  most  recent  experimental  work  concerning  the  triple  point 
of  carbon^"^  has  been  interpreted  as  indicating  a  triple-point  (melt)  tem¬ 
perature  somewhere  between  4000  and  4400°K,  it  has  been  argued  that  other 

o  19 

experimental  data  indicate  a  triple-point  temperature  near  3800  K  .  The 
corresponding  triple-point  pressures  are  argued  to  be  near  100  atm  for  the 
higher  temperature  triple  points  and  about  0.2  atm  for  the  3800°K  triple 
point. 


Resolution  of  this  very  large  (600°K)  uncertainty  in  the  melt  tem¬ 
perature  of  carbon  is  one  of  the  primary  objectives  of  the  work  presented  in 
1  2 

this  report  ’  .  A  unique  opportunity  to  do  this  is  afforded  because  the  CW 

laser  experiments  discussed  in  Section  IV  were  carried  out  using  the  same 

laser  and  essentially  the  same  apparatus  and  instrumentation  as  that  of  the 
o  19 

reported  3800  K  triple-point  experiments  .  The  importance  of  resolving 
this  uncertainty,  in  addition  to  the  desirability  of  accurate  knowledge  of 
basic  thermochemical  properties,  is  dictated  by  potentially  very  substantial 
loss  of  the  energy-dissipating  efficiency  of  carbon  materials  upon  the  onset 
of  melting^. 
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Whittaker,  A.G.,  "The  Controversial  Carbon  Solid-Liquid-Vapor  Triple  Point," 
Nature,  276(5689),  14  Dec.  1978,  pp.  695-696. 
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Noda,  T.  and  M.  Inagaki,  Bull.  Chem.  Soc.  Japan,  37,  1964,  p.  1710. 
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Fateeva,  N.S.,  L.F.  Vereshchaqin,  and  V.S.  Kolotygin,  Sov.  Phys.  Dokl.,  8, 
1964,  pp.  893-903. 
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Schoessow,  G.J.,  "Graphite  Triple  Point  and  Solidus-liquidus  Interface  Ex¬ 
perimentally  Determined  Up  to  1000  Atmospheres,"  Phys.  Rev.  Lett.,  21,  1968, 
p.  738. 
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Gokcen,  N.A.,  et  al.,  "Determination  of  Graphite/Liquid/Vapor  Triple  Point 
by  Laser  Heating,"  High  Temperature  Science,  8,  June  1976,  pp.  81-97. 


We  are  approaching  work  which  is  directed  toward  final  accurate  deter¬ 
mination  of  the  melt  temperature  in  two  phases.  The  first  phase,  discussed 
here,  is  to  place  the  melt  temperature  in  one  of  three  temperature  ranges: 
near  3800°K,  4000-4400°K,  or  greater  than  4400°K.  There  apparently  are 
no  experimental  data  which  indicate  a  triple-point  temperature  in  the  highest 
temperature  range  (see  below).  The  excellent  review  paper  by  Palmer  and 
Shelef  concludes  that  "it  seems  very  probable  that  it  (the  triple  point)  lies 
at  about  100  atm  and  4600-4800  K."  The  experimental  data  reviewed  by  Palmer 
and  Shelef,  however,  include  only  one  data  point  above  4600°K,  the  triple¬ 
point  temperature  of  4670°K  reported  by  Fateeva,  et  al.34.  Haaland3^ 
has  pointed  out  that  this  temperature  was  later  corrected  downward  to 
4040°K3^  because  of  improper  pyrometer  filter  corrections  used  in  the 
earlier  work. 


29 

Theoretical  calculations  of  Leider,  Krikorian,  and  Young  ,  based 

upon  the  carbon  species  C,  through  C7  thermochemical  data  of  Lee  and 
17  *  /  n 

Sanborn  ,  predict  a  triple-point  temperature  of  4765°K.  Theoretical  cal- 

13 

culations  based  upon  the  JANNAF  thermochemical  data  ,  which  only  considers 

carbon  species  through  C^,  predict  a  triple-point  temperature  of  about 

4900°K.  The  theoretical  calculations  do  not  consider  the  possibility  of 

22 

phase  changes  in  solid  carbon  at  high  temperature  .  Thus,  in  surmation, 
some  experimental  data  are  interpreted  to  show  carbon  melting  near  3800°K, 
most  experimental  data  seem  to  indicate  a  melt  temperature  of  4000-4400°K, 
and  recent  theoretical  calculations  predict  Tme^  >  4400°K. 


At  the  present  time,  it  is  tentatively  concluded  that  the  melt  tempera¬ 
ture  lies  in  the  4000-4400°K  range.  This  conclusion  is  based  upon  (a)  the 
weight  of  the  evidence  given  by  a  large  fraction  of  the  carbon  triple-point 
experimental  data,  and  (b)  a  preliminary  comparison  of  high-speed  motion 
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Haaland,  D.M.,  "Graphite-Liquid-Vapor  Triple  Point  Pressure  and  the  Density 
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pictures,  obtained  as  part  of  the  CW  laser  tests  described  in  Section  IV,  with 

the  high-speed  film  coverage  of  similar  tests  which  has  been  interpreted  as 

o  19 

showing  evidence  of  melting  at  3800  K  by  Whittaker,  et  al.  .  None  of  the 
most  recent  films  shows  large  "liquid"  carbon  globules,  as  seen  coming  off  the 
rapidly  spinning  sample  in  the  earlier  film  of  Whittaker  and  coworkers,  even 
though  ten  films  were  obtained  at  temperatures  above  3800°K  including  one 
for  a  measured  temperature  of  3965°K.  Thus,  the  evidence  for  melting  based 

upon  seeing  large  "liquid"  carbon  globules  is  not  repeatable.  This  could  be 

-  21 

due  to  a  number  of  experimental  difficulties  or  for  other  reasons  as  dis¬ 
cussed  below. 

It  is  acknowledged  that  other  consistently  occurring  evidence  is  given 
in  Ref.  19  to  support  the  claim  of  melting.  This  evidence  includes  spherules 
apparently  solidified  from  liquid  carbon  and  splats  on  flat  surfaces  inter¬ 
preted  as  being  caused  by  liquid  carbon  spheres  impacting  the  solid  surface. 

A  potential  explanation  of  these  is  that  they  result  from  nonequilibrium  con¬ 
densation  of  carbyne  liquid  drops  from  carbon  vapor.  It  remains  to  explain 
what  was  actually  seen  in  the  earlier  film,  interpreted  as  being  large 
"liquid"  carbon  globules,  and  not  seen  in  many  subsequent  films.  It  is  pos¬ 
sible,  based  upon  the  recent  films,  that  it  was  pockets  of  high  temperature 
carbon  vapor  which  would  have  a  viscosity  much  higher  than  that  of  the  sur¬ 
rounding  chamber  gas.  The  reason  this  behavior  has  not  been  seen  again  could 
be  related  to  an  unstable  ablating  carbon  vapor  plume  at  the  surface  of  the 
cylindrical  sample.  The  instability  may  have  been  unique  to  that  one  particu¬ 
lar  run  because  of  the  highly  irregular  and  nonuniform  ablation  of  the  pyro¬ 
lytic  graphite  evidenced  in  the  post-test  appearance  of  the  sample.  This 
interpretation  is  being  investigated  further. 

Presuming  that  our  tentative  conclusion  concerning  a  4000-4400°K  melt 

temperature  can  be  substantiated,  attention  will  then  be  focused  upon  detailed 

interpretation  of  all  the  laser  ablation  data  obtained  this  past  year.  New 

27 

pulsed  laser  tests  planned  for  the  fall  of  1980  ,  which  will  employ  an 
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auxiliary  laser  to  indicate  surface  melting  by  detection  of  the  substantial 
change  in  reflectance  upon  melting,  are  expected  to  play  a  key  role  in  the 
final  determination  of  melt  temperature.  This  technique  has  been  very  suc¬ 
cessfully  demonstrated  for  detecting  surface  melting  in  the  thermal  annealing 
?6 

of  semiconductors  .  In  addition,  we  will  be  continuing  discussions  with 
Or.  Ared  Cezairliyan  of  the  National  Bureau  of  Standards  in  Washington,  D.C., 
who  is  planning  to  carry  out  new  conventional  triple-point  measurements  for 
carbon^. 

17 - 

Cezairliyan,  A.,  Personal  Conmunication,  National  Bureau  of  Standards, 
Washington,  D.C.,  Apr.  1980. 
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